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We establish the minimum time it takes for an initial state of mean energy E and energy spread 
AE to move from its initial configuration by a predetermined amount. Distances in Hilbert space 
are estimated by the fidelity between the initial and final state. In this context, we also show that 
entanglement is necessary to achieve the ultimate evolution speed when the energy resources are 
distributed among all subsystems. 
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How fast can a quantum system evolve in time, given a 
certain amount of energy? If the system is composed of a 
number of subsystems, is entanglement a useful resource 
in speeding up the dynamical evolution? To answer the 
first of these two questions, one typically defines some 
characteristic time of the dynamics and studies its con- 
nections with the energy resources of the initial state of 
the system. Most of the previous results in this field 
0, [| HI trace back to the time-energy uncertainty rela- 
tion in the form derived by Mandelstam and Tamm [Q: 
in this way, the various lifetimes are bounded by the en- 
ergy spread AE of the system. More recently, however, 
Margolus and Levitin have pointed out that one can re- 
late the characteristic times of the system also to the 
average energy E of the initial state j5j. In particular, 
defining the lifetime of the system as the time it takes 
for it to evolve to an orthogonal state, the above results 
allows one to introduce a quantum speed limit time as the 
minimum possible lifetime for a system of average energy 
E and spread AE. In we have analyzed such a bound 
in the case of composite systems (i.e. systems composed 
by a collection of subsytems). In this context we showed 
that entanglement is a fundamental resource to reach the 
quantum speed limit when the energy resources (E or 
AE) are distributed among the subsystems. 

In this paper we extend these results by analyzing what 
happens when the quantum speed limit time is general- 
ized redefining it as the minimum time t it takes for the 
initial state g to evolve through a unitary evolution to a 
state g(t) such that the fidelity F(g, g(t)) of Q] is equal 
to a given e £ [0, 1]. Even though the scenario is more 
complex than the case e = of [||, also in this case 
it is possible to show that entanglement is necessary to 
achieve speedup of the dynamics if one wants to share 
the energy resources among the subsystems. 

In Sect. U we extend the definition of quantum speed 
limit time and derive its expression in terms of the energy 
characteristics of the initial state, first considering the 
case of pure states and then extending the an alys is to 
the more complex case of non-pure states (Sect. I A). In 
Sect. H we analyze the role that entanglement among 
subsystems plays in achieving the quantum speed limit. 
Most of the technical details of the derivations have been 



inserted in the appendixes. 



I. QUANTUM SPEED LIMIT 

Since the Hamiltonian H is the generator of the dy- 
namical evolution and defines the energy of a system, 
one expects that the energy characteristics of a state are 
closely linked to the characteristic times of its dynamics. 
In particular, we are interested on how the mean energy 
E and the spread AE relate to the maximum "speed" 
the system can sustain in moving away from its initial 
state |*). 

Take the energy basis expansion of the initial state 



|*> = E< 



(i) 



which has 
AE 



average energy E = (i$\H\ty) and spread 
tJ(9\(H - E) 2 \^). To characterizes the depar- 
ture of the system from we can use the fidelity 



p(t) = m*(t))\ 2 = 



\2 e -iE„t/H 



(2) 



where E n is the energy eigenvalue of the Hamiltonian 
H relative to \n). The quantity P(t) is the overlap be- 
tween the time evolved state \^(t)) and the initial state 
| 'I'). A measure of "speed of the dynamical evolution" 
is obtained by analyzing how fast P(t) changes in time: 
e.g. given a value e € [0, 1], how long do we have to wait 
before the state "rotates" by an amount e, i.e. before 
P(t) — e? Assuming (without loss of generality) zero 
ground-state energy eigenvalue, it is possible to prove 
that the minimum time at which this happens is bounded 
by the quantity 



%{E,AE) 



2AE 



(3) 



where a(e) and /3(e) are the functions plotted in Fig. [I]. 
Of course for e = 1 this quantity is equal to zero: in fact, 
no time has to pass to obtain P(t) = 1. On the other 
hand, since for e = we have a(e) = (3(e) = 1, Eq. (H) 



2 



reduces to the quantum speed limit time that was defined 



%{E,AE) =max( — 



7th wh 



2AE 



(4) 



which gives the minimum time it takes for a system to 
evolve to an orthogonal configuration. In the remain- 
der of the paper, however, with 'quantum speed limit 
time' we will refer to the generalized version T e (E,AE) 
Of Eq. (§). 




FIG. 1: Plot ol a(e) (continuous line) and /3(e) (dashed line) 
introduced in Eq. (g) . The insert shows the similarity between 
/3 2 (e) and a(e). 

The detailed derivation of Eq. (||) is given in App. [a|. 
Here we only give a general idea of the procedure. The 
quantity T e (E,AE) is composed of two contributions. 
The first contribution relates the speed of the dynami- 
cal evolution to the average energy E through the func- 
tion a(e) and it extends the Margolus-Levitin theorem 
||. It provides the value of %(E,AE) in what we 
will refer to as the Margolus-Levitin regime, i.e. when 
AE/E > /3(e)/a(e). The function a(e) is derived by 
introducing two functions ca^e) and a > (e) such that 
a<(e) ^ a(e) ^ a>(e). The first one is obtained by ana- 
lyzing directly the condition P(t) — e and using a class 
of inequalities that maximizes sines and cosines with lin- 
ear functions. The second one is obtained by studying 
the time evolution of a class of "fast" two level states. 
This procedure does not allow one to obtain an explicit 
analytical expression for a(e), however the two bound- 
ing functions a< and a> can be shown numerically to 
coincide giving an estimate of a(e). The second contri- 
bution to T e (E,AE) relates the speed of the dynamical 
evolution to the spread AE by means of the function 



However, for the sake of completeness, in App. |A] we 
have rederived the value of /3(e) starting directly from 
the expression (^) of P(t). 

Given E and AE, the quantum speed limit defines a 
forbidden evolution regime where the probability P(t) is 
not allowed to enter. In fact, for ^ t < Tq(E, AE), 
Eq. (||) implies 



f i (2Et\ , [2AEt\ ] 



(6) 



where a -1 and /3 _1 are the inverse functions of a(e) and 
/3(e) respectively. This regime is shown on Fig. |2[ where 
an example of P(t) trajectory is plotted. By introducing 
the Margolus-Levitin type contribution (i.e. the term de- 
pendent on E), Eq. (|6|) generalizes the previous bounds 
for P(t) jl], ||. Notice that, since the contribution de- 
pendent on AE to Eq. (^|) exhibits a null derivative in 
t = 0, it always provides a non-trivial bound to P(t). On 
the other hand, since the contribution dependent on E 
exhibits a negative slope in t = 0, it does not provide an 
achievable bound when AE ^ E. For the same reason, 



the bound ( A12 ) on the first derivative of P(t) is not 
modified by the presence of the Margolus-Levitin type 
contribution of Eq. (§: \dP(t)/dt\ is limited only by the 
energy spread AE. 




FIG. 2: Plot of the constraint given by Eq. (0), which 
shows the forbidden region where P (t) is not allowed to en- 
ter (shaded areas) . The time-energy uncertainty contribution 
to Eq. ^ defines the light grey region through /3 _1 . The 
Margolus-Levitin type contribution defines the dark grey re- 
gion through a -1 . The condition plotted here is for AE/E = 
1.73. The continuous line represents the trajectory P(t) of 
a "fast" state of Eq. ( |A7| ) with £ = 0.5. This state 

touches the boundary of the forbidden region for e = 0.30 
and* = T e=0 .3o(£,AS) =0.42 (f£). 



/3(e) = — arccos(y / e) 



(5) 



This term provides the value of T e (E, AE) in what we 
will refer to as the Heisenberg regime, i.e. when AE/E 
/3(e) /a(e). Equation (||) was previously proven in |) by 
employing the general form of the uncertainty relations. 



A. Quantum speed limit time for mixed states 

Up to now we have focused on pure states of the sys- 
tem. What happens when the system is in a mixture 
g? We will show that the notion of quantum speed limit 
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bound (U) can be extended to the density matrices in the 
sense that 7^ (E, AE) gives the lower bound to the time 
it takes for a state g with energy E and spread AE to 
evolve to a configuration g(t) such that 



F(g, g(t)) 



(7) 



where F(g, g') = {Tr [y^/gg\/g\ } is the fidelity intro- 
duced in Q. 

To prove the above statement, first of all notice that in 
the case of pure states, the fidelity reduces to the proba- 
bility P(t) of Eq. (U) and the definition (Q) reduces to the 
quantum speed limit bound given in the previous section. 
More generally, consider a generic decomposition of g, 



,Pn|</>n)(<£r. 



(8) 



where p n > 0, J2 n p n — 1 and {|0 n )} is a set of non- 
necessarily orthogonal pure states. The fidelity F has 
been shown jr| [| to satisfy the property 



F(g,g(t)) 



max l|(xlx')| 2 } 

\x),\x') 



(9) 



where \x) and \x') are purifications of g and g(t) respec- 
tively, such as the states 



n 

W) = T,vp~n\ut))\o, 



(10) 

(11) 



with {|£ n )}, {|OI being two orthonormal sets of an an- 
cillary system. Choosing \£' n ) — |£„) for all n and as- 
suming that they are all eigenstates of the ground level 
of the ancillary system, \x') becomes the time evolved 
of |x) (i-e. |x') = \x(t))) and Eq. @ implies that the 
fidelity is bounded by 



F(g,g(t))>\( X \ X (t))\ 



(12) 



Since the pure state \x) has the same energy E and en- 
ergy spread AE of g, it can rotate by a quantity e in a 
time not smaller than T e (E, AE) , as shown in the previ- 
ous section. This, along with inequality ( |l2] ) proves that 
the minimum time t for which F(g 1 g{t)) = e is bounded 
by the quantity T e (E,AE), as stated in (Q). Notice, fi- 
nally, that for e = we reobtain all the results that were 
given in [^| , since in this case the condition (Q) is equiv- 
alent to the condition Tr[^(t)^] = that was employed 
there Hi. 



where the dependence on E and AE has been dropped for 
ease of notation. Define the quantity e n = \ ((f> n \(j) n (T e )}\ 2 , 
which measures the rotation of the n-th component of the 
mixture @ at time T e . Applying quantum speed limit 
considerations to the state \(j> n ), one finds 



%(E,AE) >% n (£ n ,A£ n ) 



(14) 



where £ n and A£ n are the energy and the spread 
of the state \(j) n ). Now, among all possible purifica- 
tions of the form ( |l0| ) and ([ll]) choose one such that 
(£n|£m) = where ip nm = arg[(0 n |</> m (T £ ))]. 

From Eqs. (g) and @, it follows that 

n n 

%(E,AE) ^%(E,AE) , (16) 

where in ([TJ|) we employed the fact that a(e) and (3(e) 
are strictly decreasing functions. Combining Eqs. ( |l4] ) 
and dl^), we find that, for all n, 



T s (E,AE)>T en {£ n ,A£ n ) 



(17) 



Consider first the Margolus-Levitin regime, i.e. AE/ E ^ 
(3(e) /a(e). From Eq. (@ it follows that 



a ®2E >a{en) 2£- n 



(18) 



Since the energy of the state g is E = Yl n Pn^ni Eq. ( |is| ) 
implies 

a(e) ^ ^p n a(e n ) . (19) 

n 

Analogously, in the Heisenberg regime, i.e. when 
AE/E < (3(e) /a(e), since AE 2 = E„Pn[A£* + (E - 
£ n ) 2 }, one obtains 

^ 2 ( e -)^p/y. (20) 



The inequalities (|T^) and ( gtj ) must be satisfied if g 
reaches the quantum speed limit. Since both a(e 2 ) and 
( 3 2 (e 2 ) are stri ctly convex functions (see Eqs. ( A14[ ) and 
( A15| ) of App. A3), such conditions can be fulfilled only 
when the equalities hold: this happens if e n = e for all 
n and if the equality holds also in (|l4|) . This shows that 
the fastest states g are mixture composed by pure states 
\4>n) that all achieve the quantum speed limit bound for 
the same e at the same time. 



Mixed states that reach the bound 

Before concluding the section, let us analyze under 
which conditions a mixed state can reach the quantum 
speed limit. Assume that the state g of energy E and 
spread AE achieves the bound for a value e, i.e. 

F(g,g(T e )) = e, (13) 



II. ENTANGLED DYNAMICS 

In a preceding paper we analyzed the role of en- 
tanglement in achieving the quantum speed limit bound 
for composite systems. We found that quantum cor- 
relations among the subsystem allow the state of the sys- 
tem to evolve to an orthogonal configuration faster if the 
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energy resources are not devoted to a single subsystem. 
Here we analyze the generalized bound (||) and show that 
the same result holds even when we do not require the 
initial and final states to be orthogonal. Quantum corre- 
lations among subsystems allow the state of the system 
to rotate in Hubert space faster if the energy resources 
are not devoted to a single subsystem. 

In the following we consider the case of a composite 
system with M independent components. Its Hamilto- 
nian is given by H — ^ fc H k , where H k is the free Hamil- 
tonian of the fc-th subsystem. Since the Hamiltonian H 
is assumed to have zero ground state, we will redefine all 
the H k to have zero ground states without loss of gen- 
erality. We will start analyzing the case of pure states, 
postponing the case of non-pure states. 



A. Pure states 

Consider a composite system of M non-interacting 
parts in the initial pure separable state 



l*> = l^i)i ■■■\iPm)m , 
which has energy and energy spread 

E 



AE 



fc 



1/2 



(21) 

(22) 
(23) 



where E k and AE k are the energy and the spread of the 
state \ipk)k °f the fc-th subsystem ]ict| . The state 
reaches the quantum speed limit if, for some value of e, 
the following identity applies, 



P{%{E,AE))=e 



(24) 



where P(t) is the probability (||) of the state |^) and 
%(E,AE) is the quantum speed limit time of Eq. (||). 
For a separable state, the quantity P(t) is given by 



P(t) = P 1 (t)---P M (t) , 



(25) 



where Pk(t) = |/c(V , fe|V'fe(^))fc| 2 is the overlap of the state 
of the fc-th subsystem at time t with its initial value. 
Defining e fc = P k {%{E,AE)) and using Eq. (§|), the 
condition (E4|) can be rewritten as 



e = ei • • • e M ■ (26) 

Moreover, applying the quantum speed limit relation (||) 
to the fc-th subsystem, one finds that 



%(E,AE)^T ek (E k ,AE k ) 



(27) 



Consider first the Margolus-Levitin regime, i.e. AE/ E ^ 
/3(e)/a(e). In this case, Eq. ( p7j ) and the definitions of 
Z(E,AE) a,nd% k (E k ,AE k ) imply 



a{e) 2E >a{ek) 2E k 



(28) 



and hence, using the expression (|22j) of the total energy 
of the system, 



M 



(29) 



k=i 



Analogously, in the Heisenberg regime Eq. (|27| ) implies 
and hence 



M 



(31) 



fe=i 



A necessary condition for the separable state \^f) to reach 
the quantum speed limit is that there exists a set of e k 
that satisfy at least one of the two inequalities ( p9| ) or 
( pT| ) under the constraint (p^). Accordi ng t o the stric t 



subadditivity of a(e) and /3 (e) (see Eqs. (A16) and (A.17) 



of appendix |A 3[ ) , the relations ( p9| ) and ( |31| ) can be sat- 
isfied only when the equality holds: this happens if there 
exists a fc (say fc') such that ey = e and e k = 1 for fc =/= k' . 
Such a solution corresponds to the case in which all the 
energy resources are devoted to the fc'-th subsystem. In 
fact, for fc = fc', the relations (Eq) and (BQ) imply 



Ey > <^±E = E 
a(e) 

AE k , > ^)-AE = AE 
/3(e) 



(32) 
(33) 



where Eq. (p2[ ) holds in the Margolus-Levitin regime, 
while Eq. ( |33[ ) holds in the Heisenberg regime. Equa- 
tion (H||) and the form ( p2|) o f E require that E k > = E 
and E k = for fc ^ k' J10[. Since H k have all zero 
ground state energy, this also implies that AE k > — AE 
and AE k = for fc ^ fc'. On the other hand, Eq. @ 
and the form (j||) of AE require that AE k > = AE and 
AE k = for fc ^ fc'. 

In conclusion, the only states of the form ( |2l] ) that 
can reach the quantum speed limit (^) for some value of 
e are those in which all the energy spread AE is carried 
by the single subsystem fc'. The other subsystems are in 
eigenstates of their Hamiltonians H k . Moreover, if the 
system is in the Margolus-Levitin regime, then fc' carries 
also all the mean energy E, the other subsystems being in 
their ground states. From the dynamical point of view, 
this means that fc' is the only subsystem that rotates in 
the Hilbert space, while all the others do not evolve. 



A simple example: separable vs. entangled state 

The gap between entangled states and non entangled 
ones is particularly evident is the case in which the en- 
ergy resources are homogeneously distributed among all 
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subsystems, i.e. when E k = E/M and AE k = AE/s/M 
for all k. For the sake of simplicity we analyze an ex- 
ample in which all the subsystems are in the same state 
\^k)k = \ip)k- In this case the minimum time t e for which 
the global state |\P S ) = |i/>)i • ■ ■ \iP)m rotates by a quan- 
tity e is given by the minimum time it takes for each 
subsystem to rotate by a quantity e x ' M . Applying the 
quantum speed limit (|3J) to each subsystem, one obtains 
that t e > T e i/ M (E/M, AE/y/M). The ratio R(e) be- 
tween t e and T e (E, AE), i.e. 



R(e) > min \M 



a(e) • 



/3(eV M ) 



(34) 



shows how much slower the separable state \*ff s ) is in 
comparison with the maximum speed allowed for a sys- 
tem of the same energetic resources of \^ s ) (see Fig. g). 
From the subadditivity properties ( A16Q and (A.17), it 
follows that R(e) ^ 1. In particular, for e = 0, R(e) is 
always greater than yM, as discussed in Jof . 




B. Entangled dynamics for mixed states 

In this section we generalize the results of the previous 
section to mixed states. 

The most general separable state of M subsystems has 
the form 



Q 



£p»l* (n) )<* (n) 



(36) 



where p n > 0, J2 n Pn = 1) and 

i* ( " ) ) = i# ) )i--i^V 

with l^ l )k a state of the fc-th subsystem. As discussed 



(37) 



in Sect. I A, T e (E, AE) is the minimum time it takes for 
the state g to reach a configuration g(t) with fidelity e. 
Moreover, we already know that the state g will reach the 
quantum speed limit bound for mixed states of Eq. (Q) 
only if all the states | rotate by an amount e in 
the time T e (£ n , A£ n ) , given £ n and A£ n the energy and 
energy spread of Since |\I>( n )) is a separable pure 

state, from the previous section it follows that this is pos- 
sible only if there exist a subsystem (say the k n -th) that 
possesses all the energy resources. This means that the 
only separable states g that reach the bound (]?]) are those 
for which for any statistical realization n of the mixture 
d36| ) , a single subsystem evolves to an orthogonal configu- 
ration at its own maximum speed limit time (which coin- 
cides with T(E, AE) of the whole system). All the other 
subsystems do not evolve. Since the above derivation ap- 
plies for any expansion p n , one can say that in each ex- 
perimental run only one of the subsystems evolves. This 
is essentially the same result that was obtained in |6) , al- 
though here we considered the more general case of e ^ 0. 



FIG. 3: Plot of the lower bound of R{e) from Eq. (]3J). It 
shows that homogeneous separable states cannot achieve the 
quantum speed limit bound (given by the dashed line). Here 
M = 5. 

Consider now the following family of entangled states 



l*c) = V^-e |0)i---|0)m+£ \Eo)i---\Eo)m , (35) 

where £ G [0, 1] and |0) and \Eq) are eigenstates of 
the Hamiltonian of energy and E > respectively. 
The state |$^) represents an homogeneous configura- 
tion where each subsystem has energy S, 2 Eq and spread 
£ 2 ^/l — £?Eq. However, unlike the separable state \*5> s ) 
discussed before, for a suitable choice of the parameter 
£, achieves the quantum speed limit bound, as can 



be show n by comparison with the state |f^) of Eq. ( A7) 
in App. lAl 



Proving that homogeneous separable states cannot ex- 
hibit speedup while at least one homogeneous entangled 
case that exhibits speedup exists, we have shown that 
entanglement is a fundamental resource in this context. 



III. CONCLUSION 

In this paper we have generalized the definition of 
quantum speed limit time [g[ [| to take into account the 
case in which the system does not evolve to an orthog- 
onal state. We have used the fidelity F(g, g(t)) of as 
a measure of the "distance" between the initial and the 
final states. In this context we have analyzed the role of 
quantum correlations among subsystems in a composite 
system. As a result we have shown that entanglement 
is a fundamental resource to achieve a speedup in the 
dynamical evolution in a composite system. In fact, the 
only separable states that can achieve the quantum speed 
limit bound are those where only one subsystem at a time 
is evolving, while the others are stationary. 



APPENDIX A 



In Sects. A 1 and A 2 we derive the form of the func- 



tions a(e) and (3(e) respectively. In Sect. A3 we study 



these two functions giving some mathematical properties 
that are used in the paper. 
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1. Derivation of a(e). 

In order to determine a(e) we will: i) give a lower 
bound for it; ii) give an upper bound for it; Hi) show 
numerically that these two bounds coincide, thus provid- 
ing an estimation of a(e). 

i) A lower bound for a(e) can be obtained by observing 
that if P{t) = e, then (*]*(*)) = y/e e l6 i.e. 

^ \c n \ 2 cos = \/ecos9 



E n t 



|c„| 2 sm — 



'esmf 



(Al) 



with 9 £ [0,2tv]. Consider now the following class of 
inequalities for q ^ 0: 



cos x + q sin x 5* 1 



(A2) 



which is valid for x and where a is a function of q 
defined implicitly by the set of equations 



sin y 



y+ y/y 2 (l + g 2 )+g 2 
1 + y 2 

a(l-qy)+q 



(A3) 



f or y £ [it — arctan(l/g), 7r + arctan(g)]. The inequality 
( [A2] ) is obtained by bounding the term on the left with 
the linear function that is tangent to it and is equal to 
1 for x = 0, as shown in Fig. 0. Since we assumed zero 




FIG. 4: Plot of the inequality (|Ag) for q = f . 
a ~ .64. 



In this case 



ground state energy all the energy lev els a re positive and 
we can replace x with E n t/h in Eq. (A2). Summing on 
n and employing Eq. (|Al|), we obtain the inequality 



yfe (cos 9 ~ q sin 6>) ^ 1 — a 



(A4) 



From the definition of a(e) introduced in Eq. ([|), this 
implies 

2 



a(e) > [1 - Ve (cos 9 - q sin 9)} - 



(A5) 



Since, for a given value of 0, Eq. (A5) must be valid for 
all 9^0, then the following lower bound for a(e) can be 
obtained 



a(e) > a<(e) 



mm < max 

9 



| [1 - y/e (cos 9 - q sin 6»)] 1 1 



(A6) 



ii^ To provide an upper bound for a(e) consider the 
following family of two-level states, 



(A7) 



where ^ £ [0,1], and |0) and |i?o) are eigenstates of the 
Hamiltonian of energy and Eq respectively. The state 
|f2j) has average energy E — £ 2 Eq and energy spread 
AE = £yl — £, 2 Eq. Solving the dynamical evolution of 
the state one can show that the first time t for which 
P(t) = e is given by 



Et 



^ arccos 



e-l + 2£ 2 (l-£ 2 



(A8) 



Minimizing over £ the right term of Eq. ( |A8| ) gives the 
following upper bound for a(e), i.e. 



a(e) a>(e) 



z arccos 



e- l + 2z(l 



2z(l - z) 

where the z is a function of e defined implicitly by 
~e-l + 2z(l- z) 



(A9) 



arccos 



2z(l - z) 



l-2z 



1 



l + 4s(l 



(A10) 



Hi) The obvious difficulty in deriving the explicit 
analitic fo rm of the bounds a<(e) and ot > (e) defined in 
(|A6| ) and ( |AS| ), can be overcome by performing a numer- 
ical study of these two conditions. We will show that 
a<(e) = a>(e), thus giving an estimate of a(e). 

In order to numerically estimate a> (e) one has to solve 
Eq. ( A10| ). Using a bisection algorithm, it is simple to 
get a machine-precision accurate solution very rapidly for 
all values of e. On the other hand, the estimate of a<(e) 
requires greater care, since two different parameters — q 
and 9 of Eq. ( |A6| ) — are present in its definition. For 
each value of e i t is necessary to calculate the term on 
the right of Eq. ( A5) on a bidimensional grid of values of 
q and 9 and find for each 9 the maximum on q. The value 
of a < (e) is calculated by choosing the minimum among 
these maxima. Of course this procedure is biased, since 
it depends on the grid spacings 5q and 89. In order to 
remove the bias in the calculation result, one can repeat 
the whole procedure for different values of the grid spac- 
ing and then extrapolate the result for the spacings Sq 
and 69 tending to zero. We have used a least square lin- 
ear interpolation, where x 2 minimization allows also to 
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FIG. 5: Comparison of a<(e) and ct>(e) for some random 
values of e. The plot shows a > (e) — «< (e) with the error bars 
denoting the 89 — > extrapolation error of «<(e) (see text). 
Since all the values are compatible with zero, we can conclude 
that a<(e) = a>(e) = a(e). 



extrapolation error has been found negligible in all cases, 
meaning that a linear extrapolation is well suited. To 
reduce aliasing problems, instead of using an equispaced 
grid, it is preferable to adopt a random grid of values of q 
and 9 uniformly distributed so that the average distance 
between distinct values is 5q and 89 respectively. 

The extrapolated value of a<(e) with its error bars is 
compared with the calculated value of a > (e) in Fig. [|. 
Machine-precision accuracy is rapidly attainable in the 
calculation of a> (e) and we have considered it as unaf- 
fected by error. Since the values of a>(e) and a<(e) are 
compatible for arbitrary values of e we can conclude that 
the two functions coincide and are thus both equal to 
a(e). This allows to give the numerical estimations of 
this function that have been used throughout the paper. 
Notice, however, that a(e) is roughly approximated (up 
to a few percent error) by the function /3 2 (e) as can be 
seen from Fig. [j] 



recover an "error bar" that measures how well the linear 
interpolation works for each value of e. The error bar has 
no statistical meaning: it simply gives an idea of how well 
the linear extrapolation for 59 ^ works for the value 
of e under consideration. It can be used also to give a 
"confidence interval" for the obtained result. The 5q — > 



J 



2. Derivation of /3(e). 

The function /3(e) can be derived starting from Eq. 
by the following chain of relations: 



dt 



P(t) 



^|c„| 2 |c m | 2 (£„-£) sin 



J2\ c n\ 2 \c m \ 2 (E n ~E)e 



£ \c n \ 2 (E n E) (j2 \c m \ 2 e-^- E -W h - P{t) 

n \ m j 

i 



i{E n -E m )t/h 



(All) 



where the last identity has been obtained adding a zero 
term to the sum on n . Applying the Cauchy-Schwarz 
inequality to Eq. ( All ), we find 



dt 



P(t) 



(A12) 



which for t ^ irh/(2AE) implies f| i 



P(t) > cos 2 (^-t 



(A13) 



This means that the smallest time t for which P(t) = e is 
bounded by the quantity (3(e)-Kh/ (2A E) wi th (3(e) de- 
fined in (g). Notice that the bound ( A13 ) is achiev- 
able since, for example, the state |^j =1 /^/2) of Eq. (A7) 
reaches it. 



3. Mathematical properties of a(e) and /3(e) 

f 



Both a(e) and /3(e) are strictly decreasing functions 
(see Fig. [j]). Moreover they satisfy the following con- 
straints: 



a) The functions a(e 2 ) and /3 2 (e 2 ) are strictly convex, 
i.e. for e„ <E [0, 1], 



iV 



N 



a((X^« e ™) ) < &««(e 2 ), (AH) 

n—1 n—1 

N 2 N 



where p n > and ^2 n p n = 1- The identity in 
( A14| ) and (A16) holds only if e„ = e n t for all n 
and n! . 



b) The functions a(e) and /3 2 (e) are strictly subaddi- 
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tive, i.e. for ek € [0, 1], 



FIG. 6: Convexity condition (A14) for a(e 2 ) in the case 



N 

a(Re k 

k=l 
N 

k=l 



N 



*;=i 

N 

E 

fe=i 



/3 2 (£fe) 



(A16) 
(A17) 



with the identity holding only when there exists a 
k (say k') such that ek — 1 for all k ^ k' . 

To prove these properties, one can discuss the case of 
N = 2 and then extend it by induction to arbitrary N. 
When referred to /3(e), both propertires can be analiti- 
cally proved using its definition (j^). For a(e) we must 
instead resort to numerical verification (e.g. see Fig. ^ 
and 0). 




N = 2: A a = a(e 2 ) cos 2 <f> + a(ei) sin 2 <j> - a((ei cos 2 4> + 
62 sin 2 (j>) 2 ) > 0. In this plot e 2 = 0.7. Notice that A Q is 



null only for ei = 62 and for 



0,7T. 




FIG. 7: Subadditivity condition (A16) for a(e) in the case 
N = 2: A 6 = a(ei) + a(e 3 ) - a(6T^) ^ 0, for ei,e 2 > 0. 
Notice that Af, is null only for ei = I or 62 = f. 
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